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Abstract 

We study the quantum fermions+gravity system, that is, the gravitational coun- 
terpart of QED. We start from the standard Einstein- Weyl theory, reformulated in 
terms of Ashtekar variables; and we construct its non-perturbative quantum theory 
by extending the loop representation of general relativity. 

To this aim, we construct the fermion equivalent to the loop variables, and we define 
the quantum theory as a representation of their Poisson algebra. Not surprisingly, 
fermions can be incorporated in the loop representation simply by including open 
curves into "Loop space", as expected from lattice Yang-Mills theory. We explicitely 
O I construct the diffeomorphism and hamiltonain operators. The first can be fully solved 

^ ' as in pure gravity. The second is constructed by using a background-independent 

$_( , regularization technique. The theory retains the clean geometrical features of the 

OJj' pure quantum gravity. In particular, the hamiltonian constraint admits the same 

simple geometrical interpretation as its pure gravity counterpart: it is the operator 
that shifts curves along themselves ("shift operator"). Quite surprisingly, we believe, 
J\. • this simple action codes the full dynamics of the interacting fermion-gravity theory. 

^ I To unravel the dynamics of the theory we study the evolution of the fermion- 

gravity system in the physical-time defined by an additional coupled ("clock"-) scalar 
field. We explicitely construct the Hamiltonian operator that evolves the system in 
this physical time. We show that this Hamiltonian is finite, diffeomorphism invariant, 
and has a simple geometrical action confined to the intersections and the end points 
of the "loops". The quantum theory of fermions-|-gravity evolving in the clock time 
is finally given by the combinatorial and geometrical action of this Hamiltonian on a 
set of graphs with a finite number of end points. This geometrical action defines the 
"topological Feynman rules" of the theory. 



o 



c^ 



1 Towards Quantum Gravitational Dynamics 

1.1 Motivations 

The theoretical description of the Quantum Electromagnetic Field was de- 
veloped in the thirties, shortly after the birth of quantum theory. In spite 
of some immediate successes as the justification of the existence of photons, 
it is not untill the construction of the full quantum theory of the fermions + 
electromagnetism system, nameley QED, that the full power of this theory be- 
came clear. Recently, a certain number of advances toward the construction 
of Quantum General Relativity (GR) have been achieved ||l|, g, || ^, |§] (for an 
introduction see [g[ 0j §)• However, we suspect that this theory too will express 
its full physical value only when a fully interacting matter -I- gravity theory is 
constructed, and, in particular, when the realistic fermions-|-gravity theory is 
constructed. By analogy with QED, we shall denote the quantum theory of 
fermions + gravity as Quantum Gravitational Dynamics, or QGD. 

There is a number of reasons for suspecting that matter couplings are needed 
for clarifying the Quantum Gravity puzzle. The first of these, is that it is very 
difficult to write fully gauge invariant quantities on the phase space of General 
Relativity alone, due to diffeomorphism invariance [nl M. Equivalently, it is 
extremely difficult to imagine well-defined experiments to be performed on the 
gravitational field alone. Thus, in the pure gravity case we are in the funny 
situation of constructing a theory, but not being sure of what precisely should 
we ask to the theory - a situation quite familiar, we believe, to anybody working 
in non-perturbative quantum gravity. On the other side, diffeomorphism invari- 
ant quantities, as well as realistic experiments described by those quantities, 
can be constructed in a relative simple fashion in the presence of dynamically 
coupled matter [0, 11^, lll|. For instance, the Solar System, or a binary pul- 
sar emitting gravitational radiation, are examples of matter -f gravity systems 
that we understand well as far as measurements are concerned: we know ex- 
actly which meaningful! observables we can measure. If these observables could 
be measured with Planck-scale sensitivity, then these systems could be seen as 
quantum gravity laboratories. Clearly, we need a matter -I- gravity quantum 
theory in order to describe them theoretically. 

A second reason for coupling matter to quantum gravity, is given by the 
peculiarity of the non-perturbative quantum theory of gravity. In the Loop 
Representation |^, Q, the theory has a characteristic geometrical structure: 
Quantum states of gravity are classified by Knot Theory (I) , and the dynamics 
can be represented in a fully combinatorial-algebraic fashion on Knot Space B . 
These features are not accidental, but rather are consequences of diffeomorphism 
invariance; equivalently, they are related to the fact that the Loop Represen- 
tation is a genuine background independent quantum field theory. Since the 
theory relies on these geometrical structure, it is mandatory to check whether 



these structures are lost when further fields are coupled. If so, doubts could be 
cast on the value of the Loop Representation of Quantum GR. 

1.2 Introducing fermions in the loop representation 

Motivated by the considerations above, we have studied the fermions + gravity 
system, or QGD, in the Loop Representation. The choice of fermions is mo- 
tivated by realism, but also by the fact that they are very natural objects in 
the Ashtekar formalism. The study of the gravitational interaction of fermions 
in the light of General Relativity, goes back to Dirac and Sciania [|3|, and has 
more recently been investigated by Nelson and Teitelboim |14] using a canonical 
approach. The theory has been formulated in terms of Ashtekar variables in ref. 
[Em; we review this formalism below. We then construct the quantum theory 
following the lines along which the quantum theory of pure General Relativity 
was constructed. We define the natural extension of the loops observables to 
fermions (these are given by a parallel transport operator associated to an open 
curve contacted with fermions sited at the end points of the curve), study their 
Poisson algebra, and define the quantum theory as a linear representation of 
this algebra. In analogy with the pure gravity case, we show that the resulting 
representation can also be heuristically obtained from a naive Shrodinger-like 
representation by means of a (ill defined) Loop Transform. 

The loop representation of QGD turns out to be a very natural extension of 
the pure gravity case, obtained by including open curves into Loop Space. This 
is certainly not surprising, since the kinematics of the loop representation can 
be seen as the continuum version of the Wilson-Kogut construction in lattice 
Yang-Mills theory ^^, and from the work of Gambini and collaborators [|2| 
where fermions are represented by the end points of open lines of flux on the 
lattice. In the rest of the paper, we will denote both open and closed curves 
as loops , disregarding consistency with the dictionary. It is not difficult to 
solve the diffeomorphism constraint on the resulting state space. The complete 
classification of the solutions is given by a generalization of the Knot Classes of 
the pure gravity case - the new classes include graphs with an arbitrary number 
of intersections and open ends. Thus, quantum states of QGD admit the same 
kind of topological description as the states of pure Quantum GR, contrary to 
the fear that this aspect of the Loop Representation could be lost in presence of 
matter couplings. We view this as an encouraging result, though a result that 
could have been anticipated. 

On the other side, the results we obtain about the dynamics of the the- 
ory are unexpected and, we believe, rather surprising. The dynamics is given 
by the hamiltonian constraint, which includes the fermion kinetic energy and 
the fermion-gravity interaction. We construct in this paper the corresponding 
quantum operator, and its action on the loops turns out to have an extremely 
simple geometrical interpretation: The hamilonian constraint operator essen- 
tially "shifts" loops along their tangent. This same simple geometrical action 



of the hamiltonian constraint was recognized in the context of pure gravity in 
ref-pl. The surprising result here is that the very same action, extended to 
open loops, codes the kinetic fermion energy and the fermions-gravity interac- 
tion. This was first noticed by one of us in ref . p^ . 

1.3 Introducing a clock 

While suggestive, the above construction is not fully satisfactory for three rea- 
sons. First there is a divergence in the action of the hamiltonian operator which 
is difficult to control. Second, in spite of the simplicity of the action of the 
hamiltonian operator, we have not been able to solve the corresponding quan- 
tum constraint equation. Third, the presence of the fermions does not takes us 
completely out from the difficulties of constructing gauge invariant observables: 
it simplifies the task of finding three-dimensional diffeomorphism invariant quan- 
tities, but it docs not help with the problem of finding quantities that commute 
with the hamiltonain constraint. Thus, the actual physical content of the theory 
is still quite unaccessible, as it is in pure gravity. 

To face these problems, we take one further step. We combine our results 
on the fermions with the results obtained in ref.B. In that paper, the idea 
was proposed to unravel the dynamics of quantum gravity by coupling a scalar 
field, which could behave as a clock-field, following a long tradition |jll| of ideas 
of using matter for simplifying the gravitational theory analysis (see P, lOl 
O, pO, Bl^, E3|). It is shown in Ref. p] that by a suitable gauge fixing one 
can express the dynamics of gravity as intrinsic evolution with respect to the 
intrinsic time (or physical time) defined by the scalar field. This evolution is 
explicitely generated on the state space by a hamitonian operator H. Here, 
we extend this construction to fermions. Namely, we consider the generally 
covariant gravity-|-fermions-|-scalar field system, we solve with respect to the 
scalar field, so that the hamiltonian constraint disappears from the theory, and 
is replaced by a genuine diffeomorphism invariant hamiltonian that evolves both 
gravity and fermions in the scalar field clock time. We explicitely construct 
the quantum, hamiltonian operator H (as opposed to hamiltonian constraint 
operator), by making use of the regularization techniques on manifolds that 
have recently been introduced in quantum gravity. In the fermion sector, 
we recover in this way the simple action described above, namely the shift 
of the loops along themselves. However, now the resulting operator is fully 
diffeomorphism invariant, and finite . 

1.4 QGD 

The resulting QGD is then given by a set of quantum states represented by 
graphs with a finite number of intersections and open ends, and by an Hamil- 
tonian operator that acts in a simple geometrical and combinatorial fashion on 
these graphs. Matrix elements of this Hamiltonian can be interpreted as first 



order transition amplitudes between the graph states in a time dependent per- 
turbation expansion in the clock time. The explicit computation is complicated 
by the need of extracting the square root of an infinite matrix, a task we expect 
could be solved order by order. In the present paper, we only begin the explicit 
computation of matrix elements of the operator. 

The picture of Quantum Gravitational Dynamics that begins to emerge from 
this construction has a simple and perhaps appealing general structure: A 
graph with two open ends, say, represents two fermions interacting gravita- 
tionally among themselves, and with the surrounding gravitational field. With 
the machinery developed in this paper we could (at least in principle) follow the 
quantum evolution of this system in clock time. 

The paper is organized as follows. In section 2, we review the classical 
Einstein- Weyl theory in Ashtekar form, and we introduce the loop variables 
and their extension to fermions. In section 3, we define the quantum theory, 
and we discuss and solve the quantum diffeomorphism constraint. In section 
4, we give a preliminary discussion of the Hamiltonian quantum constraint and 
its surprising simplicity - the fully consistent construction needs the clock-field. 
In section 5, we recall the main ideas of the clock-field construction, we define 
the gravity-|-fermions-|-clock- field theory, we study the regularization of the 
Hamiltonian, and we obtain the rigorous form of the Hamiltonian operator. In 
section 6, we describe resulting structure of QGD, we collect a certain number 
of comments on the general construction, we discuss the lines along which the 
theory should be further developed, and we summarize our results. 

The signature of the spacetime metric is ( — \- H — h). Throughout units are 
used in which G = c = h = 1, except in the last section. 

2 Einstein- Weyl theory in Ashtekar form 

We consider the system formed by general relativity and a Weyl fermion field. 
We describe general relativity in terms of Ashtekar variables [g3|; and the 
fermion in terms of a two- component, massless spinor field ip (x), A = 1,2. 

A word should be spent concerning the coupling between the fermions and 
gravity. Matter spinor fields coupled to gravity in the Ashtekar's formalism 
have been studied in refs. ||l^ (see also |Q). Ashtekar et.al. have added a 
quartic terms to the minimally coupled fermion-gravity action in order to get 
the equations of motion of the Einstein- Weil torsion-free theory. As made clear 
in pq l, the minimally coupled action without quartic terms, yclds the Einstein- 
Cartan theory, in which fermions act as a source of (non-propagating) torsion. 
R. For simplicty, we consider here the minimally coupled theory. The system 

^In the theory with the quartic interaction one can naturally define a connection which is 
the sum of the torsion-free selfdual connection, and a matter generated term; the equations 
of motion are then structurally equal to the Einstein- Weil equations, but with the selfdual 
connection replaced by the full connection. See |p5|| for more details. 



represents for instance neutrinos in a dynamical spacetinie. We expect that 
the extension to Dirac fermions, and to the theory with the quartic interaction 
should be straight forward. 

In recent years it has become conventional to work with (classical) spinor 
fields which are Grassmann valued. We (reluctantly) decided to follow this 
practice because it simplifies the bookkeeping of the signs in promoting Poisson 
brackets to (anti-) commutators (since we do not use canonically conjugate 
variables as basics observables for the quantization). Thus, we assume the fields 
V' (x) anticommute. 

Spacetime is taken as a manifold M with the topology S x H, with E com- 
pact. We indicate spacetime indices (four dimensional as well as three dimen- 
sional) with lower case latin letters. Over and under tildes denote densities of 
weight -1-1 or —1 respectively (on M as well as on S). Since we are dealing with 
spinors, it is natural to adopt the spinor version of the Ashtekar formalism. The 
gravitational variables are then the spacetime 5*^(2,(0) soldering form ^(J°'^ (x), 
that is the spinor form of the tetrad field, and the SL{2X) Ashtekar connection 
[ p3| , which has only unprimed indices, ^A^^{x), and which is interpreted as the 
(internal-) selfdual part of the spin connection. We give the explicit relation 
with the vector quantities below. The Lagrangian of the system is [|6[ ^ ^ 

/:Ei„.tei„-Wc.yl = i' ^Y ^\^' ' ^\ a'Ki/'^ + ^ C of a\j,,i,^' ^V ^i^^ . (1) 

In order to construct the quantum theory, we set the theory in canonical 
form. This can be done in a fully covariant way by taking the space of the solu- 
tions of the equations of motion as the phase space, and defining the appropriate 
symplectic struture over it. However, we follow here the equivalent, but simpler 
and more conventional, approach, based on the 3-t-l decomposition, even if it 
obscures the natural general covariance of the hamiltonian formulation (see for 
instance ||2|, [29||). 

The torsion-free derivative operator compatible with ^cr°^ will be denoted 
V, while the derivative operator associated to the selfdual connection '^A^^ 
is denoted by ■*!?. Let t be a function with nowhere vanishing gradient, whose 
level surfaces Sf are diffeomorphic to E. We denote the (three dimensional) pull 
back to Yit of four dimensional quantities with the same notation as the four 
dimensional quantity, taking away the superscript "^", if present. Let also i" be 
a vector field with affine parameter t, i.e. f^V at = 1. Given a soldering form 
^<j\js^i such that each Ej is spacelike, take the normal to Sj as n" : n°'na = — 1. 
The induced metric on Ej is denoted Qab '■— gab + nam. Thereby every tensor 
field can be decomposed into its parts on and orthogonal to E(. In particular, 
i° gives us the Lapse and Shift: t° = 7Vn° -|- iV". Unprimed SL{2,(C) spinors, 
on M, and SU{2) spinors, on E, can be related as follows. One defines the 

hermitian metric in the spinor space, G^"^ := iy^n"^^ (= G^ "^) , where 
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A A' „ A A' 



^a n°. The primed SL{2,(E) spinors define a hermitian conjugation 



conjugate of ^a- The soldering form for SU{2) spinors 



operation on the unprimed ones through ("0^)4 •= G^ ipA'- The hermitian 
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a\B:^-iV2^a\/nB)A', (2) 

is hermitian {(T\g) — cr^Bi ^^'^ traceless a"'^^ = 0; it provides, locally, an 
isomorphism between TS and second rank, trace- free, hermitian spinors jab- 
In particular g"'' = -Tr {a°^a'') . 

The configuration variables can be taken as A^a s-^d ^p^, whereas the cor- 
respondent canonically conjugated momenta are —i^a°'A '■— —i\f^{o)a°'A 
and TTA := — «(o') (V"^)^- The action of the Einstein- Weyl system becomes: 

S'Einstoin-Wcyl = /^*/ ^^2:7?^ Z^t^A'^ "f Tr f-i\/2(T''j £tv4f, 

+ (t.4A^^)GAB + A^"K + iy^. (3) 

Where Lt is the Lie derivative along i". Here 

Gab = -i^PiV^t'a' AB^T^i.A'^B) 

Va = -iV2TT{^''Fab)-^A'Daij^ 

S = -TT{a''a'Fab)+iV2a%''nBVai^^. (4) 

These are the Gauss, vector and scalar constraints, respectively, for the Einstein- 
Weyl system as given in terms of Ashtekar variables. Finally, the theory is 
defined by the conventional reality conditions which pick the real sector of the 
phase space to which classical general relativity belongs. We will not discuss 
reality conditions in this paper. 

The symplectic structure one arrives at is 

{a'^AB{^),A,^^iy)} = -^S\x,y)S,^S^fSBf, 

{^A(x),V^^(y)} - -^^(x,y)5/. (5) 

The relation between these canonical variables and the vectorial Ashtekar vari- 
ables (see |6|), namely the densitized triad field E°-^(x) and the 5*0(3) connection 
A\{x), is given by 

al^(x) = -^ E''\x) cj.a'', AaA^'ix) = ^ Ki^) a^A^" (6) 

where i = 1,2,3 and at are the Pauli matrices. 

With this simplectic structure, the constraints can be interpreted as genera- 
tors of gauge transformations, spatial diffeomorphisms and the very dynamics. 



This is shown, together with the fact that they form a Poisson algebra, in 
[Ashtekar et.al. ] by smearing the constraints as 

Gt = - f d^xT^^GAB (7) 



D„ = I (fx 



v-Va - v^AJ'^Gab (8) 



Hn = iV2 £xNS. (9) 

j'AB ^ w°, N are arbitrary test fields. 

2.1 Classical fermion paths 

Non-perturbative quantum General Relativity is constructed in terms of the 
loop variables ^ 

T[a] := l/over2 U/[a] (10) 

T'^mis) := V2C//[/3](.)?V(/3(s)). (H) 

We indicate loops by greek letters. A loop is here a closed continuous piccewise 
analytical curve in S, a : S*! — > S; and s G [0, 1] is the parameter along the loop: 
a : s\^ a^{s) (We identify the values s and s + 1). We indicate by Uj^[a]{s) 
the parallel transport SL{2,<E) matrix of the Ashtekar connection around the 
loop a, starting from the parameter value s; that is, the path order exponential 
of the line integral of the connection around the loop. 
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s) := PexpiJ ds ^^ AaA^"}. (12) 



In order to construct the extension of the loop representation to fermions, 
we want to generalize these loop variables to the presence of the spinor field. We 
want to define objects invariant under the action of the Gauss constraint. Loop- 
like variables in theories with connections involving fermions, have been dealt 
with by Gambini and collaborators iQ for Yang-Mills, and by Kim et.al. pUj , 
for 2-1-1 gravity. For earlier related ideas, see |1^, ^. We follow here the lines 
of iQ. Let us consider piecewise differentiable continuous open curves in S. As 
we said in the introduction, we will call these open lines as loops, certainly with 
an abuse of terminology; and we denote them too by means of greek letters: 
a : (0, 1) -^ S; and a : s\^ a°'{s). We denote as ai and a/ the initial and final 
point of the loop, namely: 

a^ := a(0), 

af := a(l). (13) 



If aj = Pi, we denote the open line obtained joining a and /3 as a • /3; that is 

[a-(3]{s)-~a{2s), if s e [0, 1/2] (14) 

[a-/3](s) :=/3(2s-l), if se [1/2,1]. (15) 

We then define 

X[a] := V^(a,) U/[a] V^s(a/) (16) 

Y[a] := n^{a^) U/[a] V'B(a/). (17) 

X and Y, are parametrized by open curves. They are defined as path integral 
exponentials of the Ashtekar connection along these curves, with spinors vari- 
ables attached to the end points. They are clearly S'[/(2)— invariant. Other 
important properties of the X and Y variables are the following. 

1. They are invariant under a positive derivative monotonic reparametriza- 
tion of the open loops. 

2. X is invariant under inversion of the open loop, X[a^^] = X[a]. This 
important property follows from the fact the fermions are Grassman vari- 
ables. In fact: 

= -7/.^(a/)C/AB[«~V(a«) 

= -ij^{ai)UBA[a]i>^{af) 

= +V''(a«)C/s^HV'A(a/) 

= +i^^(a,)UA''[a]tPB{af) 
= X[a] (18) 

In the third and sixth line we have used the spinor index property £, pa = 
—^aP- In the fourth line we have used the parallel propagator property 
C^AB[a] = —Uba[c(^^] (recall that if the parallel propagator /7^^[a] is the 
identity, then [/^^[a] — cab)- In the fifth line, we have switched the two 
fermions, gaining a minus sign due to their Grassmanian character. 

3. Retracing identity. As their closed loops counterparts, the fermionic loop 
variables X and Y satisfies the retracing and spinor identities that follows 
from their being defined in terms of parallel propagators of an SU{2) 
connection. For instance, if (3f = ji = Si, then X[f3 ■ •j ■ j~^ ■ S] = X[f3 ■ S]. 

4. Spinor identity. The following notation is useful. Let a be an (oriented) 
open line. We define 



A 



a 



:= ^""(m) t/s^H. (19) 



Then, if aj = /?/, we can write 

X[a-p-^]^a'^f3A (20) 

Now, consider a point p where 4 hnes terminate, that is olj = j3f = ^f = 
Sf — p. There are three possible ways of connecting these four hnes to 
form two gauge invariant X variables: 

X[a-j-^]X[S-l3-'] = aA7^SBl3'' ^a^lS^'j^'S^'eAceDB (21) 
X[a-S-^]X[-ff3-^] = aAS'^-fsP'' ^a'^P^'l^S'^eADeBC (22) 
X[a-^-^]X[-f6-^] = aA(3^-fB6^ = a^p''-f^S''eABecD. (23) 

By using the fundamental spinor identity, which is at the root of the 
Mandelstam relations 

f^AB^CD + f^ADf^BC + <^AC<^DB = 0, (24) 

we have the fcrmion version of the spinor identity, namely 

X[a-r^]X[5-p-^]+X[a-5-'^]X[-i-p-^]+X[a-(i-^]X[-i-5-^] = 0. (25) 

Similarly, it is simple to derive the identitity that refers to the intersection 
between the open loop a-P and the closed loop 7, where af = (ii = ji = 7/: 

X[a ■ /3] T[7] = X[a • 7 • /3] + X[a ■ 7"^ • /3] (26) 

5. Fermionic (Grassmann) identities. Since Grassman variables anticom- 
mute, and since the fermion field has only two components, if we multiply 
three or more fields in the same point we obtain zero. It follows that 
we can have products of X variables with at most two coinciding hands. 
Thus, for instance, if ai — Pi — ji, then X[a]X[/3]X[7] = 0. 

6. Gauge observables algebra. Finally, the most important property of the 
X, Y, T, T"' variables is that their Poisson algebra closes. A direct compu- 
tation yelds 

{X[l3],X[a]} = 

{Y[l3],X[a]} = 6^{af,l3,)X[a-l3]+5''{l3,,ai)X[a-^-l3] 

{Y[alYm = 5^{af,p,)Y[a-[3]-5^a,,Pf)Y[P-a] (27) 

Whereas the nonvanishing brackets with the T variables are 

{r«[7](s),XH} =. *A'^[7(s),c.] Y. M^b#.7l 

{T'^[7](s),r[a]} = zA'^[7(s),a] ^ M?[a#.7l- (28) 



As in the pure gravity, one may define also higher order observables, which 
will have simple quantum operators associated. For instance we can insert 
"hands" into the X variable, and so on. In order to treat the dynamics, the 
following variable will be particularly useful. 

Y-[a](s) :- ^^(aOC^A''H(0,s)a'^B^(a(s))[/c^[a](s,l)^z?(a/), (29) 

where the notation [//i^[a](s, t) indicates the matrix of the parallel trasport 
along a from s to t. This variable is quadratic in the momenta and will play 
a role analogous to the role of the T^ variable in pure gravity, which is also 
quadratic in the momenta. Indeed, we will use it for defining the hamiltonian 
constraint. 

Finally, we may introduce a further small simplification in the notation: if 
the context is clear, we may write X[a\ as T[a\. Namely we can use the same 
notation, T[a], to indicate the the trace of the holonomy of a if a is a closed 
loop, and to indicate the parallel propagator along a sandwhiched between two 
fermion fields if a is open. 

3 QGD: kinematics 

Let us now begin the construction of the quantum theory. Following the phi- 
losophy described in refs. |^, ^ |34J, we look for a quantum representation of 
the classical loop algebra. Rather than simply writing the rappresentation, it is 
perhaps more instructive to use the Loop Transform, introduced in ref. 0] , as an 
heuristic device to help us in this task. Thus we first consider a " Schrodinger- 
like", rapresentation of the quantum Einstein- Weyl theory |^. We consider 
functionals ^[A, ^] on the configuration space, and we define the canonical co- 
ordinates operators A and ip as multiplicative operators, and the corresponding 
momentum operators as functional derivative operators 

Since we are using this construction only as an heuristic tool to find a possible 
form for the fermion loops operators, we are not particularly concerned with 
mathematical rigor here. We focus on loop states in this representation. These 
are defined as follows. Given a closed loop a, we write 

^a.[A,i>]=U[a]A^. (32) 

Given an open loop a, we write 

^4A,^]^^lj'^ia,)U[a]A''^Biaf). (33) 



10 



Given an arbitrary collection /3 of a finite number of open or closed loops 
ai, a2, ..., a„, we write 

^,3 =*ai'J'a2--*a„- (34) 

Note that we follow here the usual convention of denoting loops (open or closed) 
and multiple loops (namely collections of a finite number of loops) by means of 
the same notation, that is greek letters from the begining of the alphabet. The 
idea of the loop representation is to take the '^p states as an overcomplete basis 
of quantum states. We thus introduce the Loop Transform Q as follows. 

*[/3] := fvAVtP */3[A,7/>] *[A,V]. (35) 

For a rigorous mathematical definition of these kind of integrals see the recent 
work of Ashtekar and Isham, and Ashtekar and Lewandowsky [B6| Here ^[A, V'] 
represents a generic wave functional in the connection representation, and ^[/3] 
represents its Loop Transform. The novel features entering here come from the 
inclusion of fermions; states containing no fermions are described exactly as in 
the pure gravity case W . The transform (pq) has a well defined meaning on the 
lattice [?], where the ^^ states are the Wilson-Susskind states. In the lattice 
case it is possible to show that the transform defines a unitary transformation 
to a new basis in the Hilbert space of the theory. 



Eq.(|35[) suggests that we may look for a representation of the full Weyl- 
Einstein loop algebra on a space of functionals of multiple loops, where the 
multiple loops are sets of closed as well as open loops. Then the transform gives 
us immediately the action of the X and Y operators in the loop representation 

(x [a] *) [/3] = / VAVtP *0 [A, V'] (x [a] *) [A, i^] 



= *[/3Ua]; (36) 

where the union operaton U of set theory is well defined between the multiple 
loop P and the single open loop a. Thus we have 

X[a]*[/3] =*[/3Ua]. (37) 

Note that this is the essentially the same action as the action of T[a] in pure 
gravity. Using the notation suggested as the end of the previous section, we can 
write 

f[a]*[/3] =*[/3Ua] (38) 

for open as well as for closed a's. 

11 



As far as Y is concerned, we have 

Y[a]-i')[f3] = jVAD^j'$p[A,i'](Y[a]-i')[A,tP] 



= Y,5\a,.Pf)no^-(i]+Y.5\a,,p,)%f3-a-\ (39) 

Here X]/3 indicates the sum over all the initial points of the open loops in the 
multiple loop /3, and ^^ . indicates the sum over all the final points of the open 
loops in the multiple loop (i. We introduce the following notation. We write /3e 
to indicate any end point of the multiple loop /?. If /3e = ai, the notation (3 -e a 
indicates the multiple loops obtained by attaching a^ with P^. Thus we have 

yH*[/3]=^53(a„/3e) *[/3-ea]. (40) 

/3c 

In words, the action of the operator Y[a\ is simply to attach the open loop a 
to any open end that happens to be in the point aj. 

Next, we supplement (|3^), (^0|) with the usual quantum T- variables in the 
loop representation [nl o]. The computation of the quantum commutation re- 
lations of the entire set is then straigthforward. The result is that the set of 
operators X,Y,T,T^ provides a representation of the classical algebra (p7|), 
( p8| ) and the classical T-algebra. 

We can also naturally introduce quantum operators corresponding to higher 
order loop variables m. As their pure gravity counterparts, these have a quan- 
tum algebra that reduces to the corresponding classical Poisson algebra in the 
limit in which the Planck constant goes to zero. We write here the quantum 
operator corresponding to the y° variable defined above, since it will be used 
in the construction of the Hamiltonian 

Y''[a](s)np] = y^6\(i,,a,) [dtp%t)5\a{s),m) 

Here we have indicated by a * >i=^f/3 and a * *Jt/3) the two loops obtained by 
joining the /3e end point of j3 with the point a^, and rerooting the intersection 
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a{t) — P{s) in the two possible ways. Note the plus sign between the two terms, 
which will play an important role in what follows. The expression mV) can be 
obtained for instance from the loop transform. To determine the correct overall 
coefficient and sign, an accurate computation with the SU(2) index algebra is 
needed. Note however, that the relative plus sign between the two terms in 
parenthesis is forced by symmetry, since neither of the two can be preferred. 

3.1 Diffeomorphisms and diff-invariant states 

The classical vector constraint generates spatial diffeomorphisms when acting on 
gauge invariant objects. This is also true in the quantum theory providing the 
correct ordering of the vector constraint quantum operator is chosen. Precisely 
as in the pure gravity theory, there are several equivalent ways for reaching this 
conclusion: 



1. As suggested by Isham 34 , the vector constraint can simply be defined 
as the generator of the natural action of the diffeomorphism group on the 
space of the open and closed loops. The commutator algebra of these 
generators among themselves and with all the other operators in the the- 
ory, then, reproduces the corresponding classical Poisson algebra. This 
is a sufficient condition to insure that the classical limit of the quantum 
theory that we are constructing reproduces the classical theory we started 
from. Since the correct classical limit is the sole requirement we have 
on the theory, the quantum diffeomorphism constraint defined in this way 
represents a consistent quantization of its classical counterpart. 

2. We can use the transform, and define the loop representation vector con- 
straint opererator as the transform of the vector constraint operator in 
the representation that diagonalizes A and ip. 

3. We can express the classical vector constraint in terms of loop variables, as 
a suitable limit of a sequence of these variables. The corresponding quan- 
tum constraint is then defined as the limit of the corresponding quantum 
loop operators. 

As in pure gravity, it is not difficult to show that these different strategies yeld 
the same quantum diffeomorphism constraint, and that this is can be expressed 
as follows. For every diffeomorphism (j> G Diff[E] 

*[a] = *[(/)• a] (42) 

where [(f) ■ a](s) := 4>{a{s)). The general solution of the diffeomorphism con- 
straint is the given by the loop functionals constant along the orbits of the 
action of the Diff group on the loop space, namely, they are given by 

*[a] = *[if(a)], (43) 
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where K(a) is a generalized knot class, that is, an equivalent class under difFeo- 
morphisms of sets of graphs formed by open and closed lines. 

For every generalized knot class K, we can define a corresponding quantum 
state ^K as the characteristic function of the class. We will also use a Dirac 
notation '^[a] — (a|^), and denote the state ^^ as \K >. Thus 

{a\K) = 1 if K = K{a), 

{a\K) ~ otherwise. (44) 

Let us begin here some preliminary investigation of the structure of the ensemble 
of quantum states \K >. 

Consider a fixed class K. Let a be one of the (diff-equivalent) multiple loops 
that belongs to K. Let a be composed by c closed loops and o open loops. 
There are 2o end points a^ in a. We distinguish the end points as simple or 
doubles. An end-point ag is simple is there is no other end point a'g in a such 
that tte = ci'e- It is double otherwise. Note that there may not be "triple" end 
points, because of the Pauli principle. Let S and D be the number of simple 
and double end points. And N = S + 2D = 1o be the number of end points. 

Consider the points i in the immage of a such that one or more than one of 
the following is true: 

1. a is non injectivc in i, 

2. i is an endpoint, 

3. a is non differentiable in i; 

we denote these points i as generalised intersections, or, simply as intersections. 
We assume that the number of these intersections is finite, and we denote this 
number as /. Given an intersection i, we assume there is only a finite number 
of components of a coming out of it. We denote this number as m^, and we call 
it the order of the intersection. Intersections of order 1 are single end points. 
Intersections of order 2 are either double end points or kinks along a loop. 
Intersections or order 3 are single end points that fall over a loop. Intersections 
of order 4 are either crossings of two loops, or a double endpoint that falls over 
a loop, and so on. We call the intersections of order one free end-points. 

Consider an intersection i of order rui in a loop a. Let aj{s), with j — l...mi, 
and 0:^(0) = i, he the nii lines (components of the loop a) that come out from 

i. Let Ij, j ~ l...mi be the mi tangents of aj{s) in i. Let then ]}■ , where k is 
a positive integer, the fc— th derivative of the j—th component of the loop in i. 
For instance, ii mi — I, and i = a(0), then li ~ d 
vecalpha{s) / ds , and Zj — d^a{s)/ds'^, and so on. 

The vectors ZJ transform among themselves under reparametrization of a 
and under diffeomorphisms. Consider the space of all possible intersections of 
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order rrii. Consider to of these intersections as equivalent if they can be trans- 
formed into each other by diffeomorphisms or reparametrizations. Denote the 
space of the resulting equivalence classes as the moduli space of the intersection 
of order m^. The moduli space of an intersection of order n < 5 is discrete; not 
so, in general, for larger n. However, the moduli space is always finite dimen- 
sional. Let al' be a collection of parameters that coordinatizes the moduli 
space of the i intersection, as well as characterizing the roofings of a through 
the i intersection. 

Finally, let JCp be a discrete index that labels the braids with P (ordered) 
open hands and c closed loops. We can then (over-) characterise a quantum 
state as (see ref. [^): 

\N,I,D- aT\...aT'; IC°j^ ,J (45) 

where, we recall. A'' is the number of end-points, / is the number of intersec- 
tions, D is the number of double end points, mi.... mi are the orders of the I 
intersections, a]^^' ....aj'' are the moduli space parameters of the intersections 
and /C^ is the discrete topological class of the braid obtained by deleting 

all intesection points from the loop. At the end of the paper we will describe 
some simple generalized knot classes and make use of this notation. 



4 The hamiltonian constraint has a simple ac- 
tion 

Our last and main task is to deal with the dynamics, which is contained in the 
hamiltonian constraint. (The Hamiltonian constraint of the pure gravity Loop 
Representation was discussed in [|l|, |l2|, ^. For a comprehensive and critical 
discussion of the various approaches see ref. |3^.) We shall perform this task 
in two stages. In the present section we introduce a simple and naive non- 
regularized definition of the quantum hamiltonian constraint. This is not really 
satisfactory because it does not allow us to control the divergences of the theory. 
However, we think it is usefuU to present this non-regularized version of the 
dynamics first, because it allows one to appreciate the striking simplicity of the 
geometrical action of the Weyl-Einstein hamiltonain constraint, which otherwise 
could improperly appear as an improbable product of the technicalities of the 
regularization procedure. In the next two sections, we will transform the formal 
esxpressions we obtain here in a more solid result. 

Thus, we begin by defining the hamiltonian constraint by using the simplest 
procedure: we define it in the connection representation and we transform it to 
the Loop Representation by using the Loop Transform ||l[ . We choose the order- 
ing in which the momenta are always to the right of the configuration variables. 
Taking into account the Grassmann character of the fermionic variables, the 
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Hamiltonian constraint, smeared against a test (inverse density) scalar N{x) is 
H[N] = - f (fx NTt (Pat^^aj - i V2 Tf Vatp^ ^s 



^ ^ d'x N{=^) FaM^) ^ ^ 



-^Vai;^ix)-J-—-^ (46) 

In order to compute the Loop transform of this operator, we have to compute 
its action on the kernel of the Loop Transform, that is, on the basis loop states 
\E'q[A, ^p]. We need to compute 

H[N]-i'^[A,^]. (47) 

If a is formed by closed loops alone, then Vl/a[A, ^] is independent from ip and 
therefore the second term in ( |56| ) does not act. It follows immediately that the 
action of H[N] on the closed loop states is fully equivalent to the action of the 
pure gravity hamiltonian constraint. Let us then assume a is a single open loop. 
A straigth forward calculation gives 

H[N] *„[A,^] = 



a Jo ^ 
i^''{a,)UG[F{Q.t)F,/^{a{t))UB]H{t,l)il^"{af) 
1 
ds N{ai)5\a{s),ai) a%Q) VaV{ai)UAE[a] V'''(a/) 



+ / dsN{as)6''{a{s),af)a''{l)i}'^{a^)UDA[a\Va^^{af). (48) 
Jo 

We immediately see the difficulty in this equation: the three- dimensional delta 
functions are integrated only against two line integrals, leaving a divergent fac- 
tor. As we said, let us disregard this infinity for the moment. 

We now recall from ref.|]l| that, if we disregard divergences, the action of the 
Hamiltonian operator on a pure loop state can be written as 

H[N] *„ = S[N]^^ (49) 

where iS, denoted as the "Shift operator" is a simple operator acting on the loop 
argument , as 

S[N]^o.:= ds / dtN{a{s))a''{s)5\a{s),a{t))——--^^. (50) 

Jo Jo da'' it) 
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If a has no self-intersections, and up to a divergent factor k, the Shift operator 
becomes simply 

5[iVl*a = k lim *a,„ (51) 

^^' t->0 rL ^ ' 

where 

a^Nis) := a"(s) + tN{a{s))a''{s) (52) 

The action of the Shift operator has thus a very simple geometrical interpre- 
tation: it shifts the loop ahead along its tangent. Clearly, it sends a smooth 
closed loop into itself, while it deforms a loop with kinks or intersections. This 
simple action is one of the ways in which one can (formally) understand the well 
known result that loop states corresponding to smooth non intersecting loops 
are in the kernel of the Hamiltonian constraint. The fact that the action of the 
hamiltonian constraint on the kernel of the Loop Transform can be expressed in 
terms of an operator acting on the loop variable enables us to interpret this op- 
erator as the operator that represents the Hamiltonian constraint in the Loop 
Representation R]. Thus in pure gravity the Hamiltonian constraint can be 
simply expressed as the Shift operator. 

How does the picture change if we include the fermions ? The striking result 
that we have mentioned in the introduction is that the picture does not change 
at all with the introduction of fermions: Equations (p9[), which expresses the 
Hamiltonian as the Shift operator still holds. The shift operator is still given by 
equation (|5C|), where now we also allow the loops to be open, or in the absence 
of intersections, by equation (pl|). Indeed, by computing the action of 5 [TV], as 
defined in equation ( pO| ) on an open loop state, we get precisely the right hand 
side of equation (|48| ) . What happens is that the fermion term in the classical 
hamiltonian constraint give rise to the second and third term in equation (Esf) 
and these terms are precisely the terms that "move" the two fermions at the 
end of the loop in the correct direction ! 

Thus we have the following result: 



• 



In pure gravity the action of the hamiltonian constraint on Loop Space 
can be expressed as the action of the Shift operator (pG), which simply 
shifts non- intersecting loops along their own tangent. 

• By applying this same geometrical operator on open loop states, we have 
the action of the Weyl-Einstein hamiltonian constraint. 

This result was first obtained by one of us in rcf. ||l^ . The present paper rep- 
resents an effort to clean up this result from the difliculties due to the divergent 
factor k. 

We consider the fact that the coupled Einstein- Weyl equations can be coded 
in the extremely simple Shift operator as a striking fact. We have been deeply 
puzled by this result, and we do not see any simple way of interpreting it. We 
have not been able to find any reason for which this result could be understood in 
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terms of the classical fermion dynamics. Its sim.plicity seems to us an indication 
of something, but we have not been able to decode the indication. We will return 
to this discussion in the conclusion section. For the moment, let us study how 
we can free the result from the divergences. 

5 QCD: dynamics 

5.1 Coupling a clock field 

Following the ideas introduced in the introduction, we now modify our point of 
view, and consider a richer theory: We couple a scalar field to the Einstein- Weyl 
theory. Our aim to use the scalar field in order to define a physical internal time, 
or clock-time, in terms of which we can study physical evolution. 

We have two independent motivations for chosing this roard. First this 
procedure allows us to unravel the physics of the general covariant quantum 
theory, otherwise hidden in the frozen-time formalism, as discussed in detail 
for instance in refs.y, |lO[. Second, this is a way to overcome the divergence 
difficulties we had in the previous section. Indeed, we have learned from the 
experience in pure Quantum GR that non-diffeomorphism-invariant operators 
tend to be ill-defined in a generally covariant quantum field theory, while all 
the diffeomorphism invariant operators that we have been able to construct 
have good finiteness properties 0. Thus, we wish to replace the Hamiltonian 
constraint with some diffeomorphism invariant operator [the Hamiltonian con- 
straint, being a scalar (density) is diffeomorphism covariant, not diffeomorphism 
invariant] . As shown in rcf . [gj , the coupling of a scalar clock field and the re- 
placement of the hamiltonian constraint with an hamiltonian is a way to acheive 
this result. The hamiltonian that generates the evolution in the clock time is a 
diffeomorphism invariant quantity and replaces the hamiltonian constraint. 

We refer to g for the details of the scalar field construction and the gauge 
fixing that allows to define the hamiltonian. We simply recall here the main 
idea, so that this paper could be idependently read. Physical quantum states 
are represented (say in the connection representation) as functionals ^[^,1/'] of 
the spacial fields A(£),-0(a^), satisfying the Wheler-DeWitt equation. As well 
known, the time coordinate t disappears from this frozen time formalism. In 
theory, the disappearence of the coordinate time is not a problem, since the 
observables of the theory must be 4-dimensional general covariant anyway , and 
thus, in particular, must be independent from t anyway. Examples of these 
4-dimensional general covariant observables are given by the invariant distance 
dp of the solar system planets from the Earth, seen as a function dependent 
from, say, the invariant distance d of the Earth from the sun, or as a function 
of the proper time (on the Earth) lapsed from some initial fixed event. These 
observables are manifestly coodinate invariant, and thus independent from t. 
In practice, however, it is notoriously too difficult to write the dynamical equa- 
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tions of GR directly in terms of coordinate invariant quantities: in classical 
GR we almost always work with coordinate dependent quantities, and extract 
coordinate invariant predictions only after the dynamics has been fully worked 
out in a particular gauge. For instance, we study the motion of the planets, 
including perhaps emitted gravitational radiation, in some arbitrary coordinate 
system (with some arbitrary coordinate time t: dp{t) and d{t)), and only when 
the dynamics has been solved we compute coordinate invariant quantities dp{d), 
which can be compared with astronomical observations. The quantum frozen 
time formalism, however, does not allow us to work with quantities dependent 
on the fictitious arbitrary coordinate time i, and thus makes the dynamical 
analysis particularly cumbersome: what one should do is to view the physical 
states ^[A,-)/)] as coding the quantum evolution of any variable in terms of any 
other variable. In general this is not easy. 

The problem can be simplified by studying a version of the theory in which 
there is a simple quantity to be taken as the independent variable; that is, in 
which we know from scrath which variable we want to use as the " clock variable" . 
We thus introduce a scalar field T{x, t), and we decide to study the evolution of 
the gravitational and fermions degrees of freedom, as they evolve in the value 
of T. If A{x,t),tp{x,t),T{x,t) is a solution of the field equations, we extract 
information invariant under a time coordinate transformation by solving t with 
respect to T, and substituting the resulting t{x, T) in A and ijj: We get the two 
functions 

A{x,T) - A(x,t(f,T)), 

V'(x,T) = i,(x,t(x,T)), (53) 

which are invariant under coordinate time reparametrization. Equivaletly, we 
choose a (physical) coordinate system in which r(x, t) = t. Of course, this 
cannot be done in general. It can only be done for certain solutions of the field 
equations, and for certain regions of spacetime. Thus, we are now beginning 
to use a formalism that holds only within a certain regime, and breaks down 
in approaching the limits of this regime - for instance, when the derivative 4^- 
vanishes. 

In the quantum theory, the frozen time formalism is defined in terms of 
the functionals ^[A, ^, T] of the spacial fields A(x), ^(x), r(x). We can inter- 
pret these states as giving the amplitude for a A{x),ip{x) configuration at the 
given configuration T{x) of the clock field , and to interpret the Wheeler-DeWitt 
equation as an evolution equation in the multifingered time T{x). We can then 
further fix the multifingered T{x) time as a constant function T{x) = T, (we 
keep the same letter T to indicate also the real number T , beside the function 
T(f)), and restrict -^[A.iJj.T] [functional of the field T(x)) to *[A,'(/'](r) = 
\I/[^, ■)/), T]|t(S)=t {function of the real numeber T), without any loss of informa- 
tion. The state ^[^,i/'](T) expresses the quantum amplitude for the evolution 
in T of the fields A{x, T), V'(x, T) defined above. Moreover, if *[A, V', T] satsfies 
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the Wheeler-DeWitt constraint, it is shown in [pj that ^[A, ?/'](T) satisfies the 
equation 

1%^ nA,i>W) = H *[v4>](T), (54) 

where the operator H will be defined in a moment. We can view equation (B3) 
as a genuine Schrodinger equation which evolves in the time T . 

More precisely, we can express the Einstein- Weyl-scalar-field theory in a 
canonical (gauge-fixed) form, in terms of the configuration variables A^^^ and 
■!/; , the gauge and diffeomorphisms constraints 

Gab = -iVT.Vb'a-^j^g+T^i^Ai'B) 

Va = -iV2Tr {a'' Fab) -HAVai^^ (55) 

and the hamiltonian (as opposed to hamiltonian constraint) 

H^ I (fx y'-Tr {d'^a^Fab) + iV2a''/ ^B^a^^. (56) 

It is shown in B that the solutions of this theory, A{x, T), i^{x, T), are related to 
the solutions of the full Einstein- Weyl-scalar-field theory via equation (p3|) . We 
can define the kinematics of the quantum theory and treat the diffeomorphism 
constraint exactly as in the previous sections, but now we do not have a quantum 
hamiltonian constraint, but rather a Schodinger equation (p3) and a quantum 
hamiltonian H, which is the quantum operator corresponding to the classical 
variable H given in (pq). 

The rest of this section is devoted to the construction of the quantum opera- 
tor H. This time we will not be content with formal manipulations of divergent 
expressions, and we require a somehow higher level of rigor. 

5.2 Definition of the regularized hamiltonain 

In this section we derive the key techincal result of the paper, namely the explicit 
regularized action of the Hamiltonian on the knot classes; the section is quite 
long and technical - the final result is summarized in the following section 5.3. 
Following ref. |^ , our first step is to definine a regularised classical hamiltonian 
in terms of loop variables. We introduce a fictitious background flat metric and 
a preferred set of coordinates in which this metric is euclidean; in the following 
everything is written in those coordinates and we use the euclidean metric to 
define norms of vectors and raise and lower vector indices. It will be our task, 
later, to show that the operator we define is independent from the regularization 
background metric introduced here. We write 

H= lim HLSAr. (57) 

L^O (5^0 
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HlSAt = X! -^ V ~*"Ei'ifstein / - C'^cTy* J . (58) 

Here we have partitioned three-dimensional space into cubes of sides L, labeUed 
by the index /. The quantity C^j'^^gJ-^jj^ j is the regularized form of the pure 
gravity hamiltonian constraint (integrated over the /— th cube), and has been 
defined in |5|. 

With the purpose of defining the fermion term C^J"'j j , we begin by defining 
the open loop 7^^, where r is a regularization parameter (to be taken small), 
:r is a point in space and y is a vector in the tangent space around x. Since 
we have a background flat metric we can write expressions as a; + sy , for any 
real s, with obvious meaning. The open loop 7I - is defined as the (uniformly 
parametrized) straight line (in the background metric) that starts at x in the y 
direction and is long r. That is 

{llyTi^) = x^ iilyT^s) = ^y\ (59) 

Note that 

(7s.,T(l) = ^'' + ^2/" (60) 

and 

{ilyT{\y\lr)^x'^ + y''. (61) 

By making use of this loop, we define the regularised fermion hamiltonian con- 
straint by 

-iL,r,5 ^ / j3 /^T-S 



Cwcyl / — JJ '^ ^ ^Woyl(^) ' (62) 

C^eyi(^) = CrsJd'yO{6-\y\)^^Y^[^Z^]{\y\/T). (63) 

crs = -T^- (64) 

Here 0{x) is the conventional step function, that is the characteristic function 
of R'^. Note the three different roles of the three regularization parameters 
L, T and S: the parameter L fixes the size of the boxes. As we will see, the 
introduction of these boxes will allow us to deal with the square root. For 
every space point x, the fermion term of the Hamiltonian constraint CwGyi(a^) 
is approximated by means of the loop variable Y"" corresponding to the " small" 
loop 7^-, that starts at x. The parameter r gives the length of the "small" 
loop. The direction of this loop is integrated over {(Py angular integration). 
¥"■ has a special point where the hand is inserted. This point is chosen to 
be 7j^(|y|/T) — x"" + y". Thus, also the position of this point (the hand) is 
integrated over {(fiy radial integration). This smearing of the position of the 
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hand is what determines the point spUt of the functional derivatives operators in 
the quantum theory. Note that the cPy integral is restricted (by the 9 function) 
to a ball of size S around x. Thus, the parameter S gives the point splitting 
separation between the initial point of the loop and its hand. Note that in order 
this definition to make sense we must have t > 6. 

A straigthforward expansion in r and S shows that equation ( p7| ) is satisfied, 
namely that the quantity that we have defined represents a genuine regular- 
ization of the Hamiltonian. For completness, let us sketch this expansion. We 
begin by writing the regularized expression explicitely, by using the expression 
for y°, given in eq.(p9|) 

Cwcyl(^) = CrS Jd'y 9(6 - \y\) ^ TT^hk.WA'' [l^Sy] (0, lylM 

a'^B^'iYsyilylMPc'' [Ysy] {\y\/r, i)Vi?(7.V)- (65) 

By using the explicit definition of the loop 7I-., this becomes 

Cwcyi(^) = crsjd'y 9(5 - \y\) ^ 7r^(f)[/^^ [Y,y] (0, \y\/T) 

a'^B^(x + y)f/c^ [7^,^ {\y\/T,l)i;D{x + Ty/\y\). (66) 

We may now expand everything in powers of S and r around the point x. Before 
doing so, however, let us note what follows. Since we have a {(fiy) integration 
over a volume of order S^, and we divide by S^ (in Crs) before taking the limit 
6 —^ 0, only the term in the integrand of zero order in S may survive the limit. 
We also divide by r and take the limit r -^ 0, thus only terms of first order in 
r survive in the limit. This means that to the relevant order we may replace 
quantities in x + y by quantities in x (recall the d^y integration is over a sphere 
of radius 6), and quantities in a; + Ty/\y\, which is a distance r from x by the 
first two terms in their Taylor expansion around x. By doing so, the first of 
the two parallel propagators is just replaced by the identity, while the second 
can be replaced with the entire parallel propagator along the small loop. We 
therefore have, up to terms of order S or r, 



C^fyi(^) = CrS d^y9{6-\y\)^ TT^'ixjS 



A 



\y\ 

a'^B^ix) {6c'' + ry^A^ci^)) ii'Dix) + Ty^dc^oix)) = 
= Crs [d'y 9{S ~ \y\) ^ n^{x) 

J \y\ 

a%^(f) [i^c{S) + Ty- {A^c{^)i^D[x) + dci^c{x))] (67) 

The first term in the square brackets vanishes because J d^yy"" — 0, and the 
second term is the covariant derivative. Thus 



^WoyU^) — ^■rS''' 



d'y 9{S - \y\) y^ 
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TT''{x)a''A^{x)D,^Jc{x). (68) 



The integration is now immediate 

Jd^y 9(5 ~ \y\) y^ ^ 3 S^^ ^nS^ (69) 

Restoring the exphcit form of c-rs wc thus obtain 

C^cyi = TT^a-A^^D.^c + 0{S) + 0{r), (70) 

namely the fermion term of the hamihonian constraint, as we wanted. 

We have shown that the quantity ( |63| ) provides indeed a regularized form of 
the fermion hamiltonian constraint. Wc now come to the quantum theory and 
define the corresponding regularized quantum operator 



HlSAt — 7 . ^ V "^Einstein I " ^wJ^l I ■ C^^) 



I 

~iL,T,5 ^ I i3„, At.S 



c^:^i - jijfyC'^Uv)^ (72) 

C'^eyi(y) = ^rs j d\ 9(6 - \y\) |j Y^ [jZy] (lyl/r) . (73) 

The operator ¥"■ is defined in equation (El]). Note that in this regularized 
operator the two hands of the operator |1[ do not overlap: they are point split. 
Let us study the action of the operator that we have defined: in spite of it 
apparent complexity, this action will turn out to be relatively simple. Plugging 
the explicit definition of y (El|) in d73) gives 

C^Jyf (^) ^H = CrS fd?y 6(6 _ |y|) ^ ^ S^X, ae) 

/ ds a'^is) S^Y,y{\y\M, a{s)) ^ * [« * *«,,,7:,](7.4) 

We now use the explicit form of the small loop, and we keep only terms small 
in r and 5. The first 5^ function in the last equation forces the loop a to have 
an end-point ae in x. The integration over the small ball (size 5) around this 
point, and the second 5^ function, pick up a second point a{s) — x + y in a, 
close to ae- If ae is a free end point, we can write this second point (up to the 
relevant order) as x + \y\le, where le was defined in the previous section as the 
tangent of a in the end-point ae- In general, however, ae needs not be a free 
end-point; in the general case there will be several components of a originating 
from Ug. Thus le takes a finite number of values, and the radial d'^y integration 
together with radial part of the second delta function pick up all these values; 
note that they turn out to be proportional to y. The radial d^y integration is 
then straightforward, and we obtain 

^^:yf (-) *H = JJ^i: ^^("e, X) E E * [[" * -'^■^%i,+Sl)\ • (75) 

Oe ;-; q=± 
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We introduce the notation 

c^:;^i^) *H -jjsj2 '^'("-^) ^e' *N (76) 

^I' *H -T.T.'^ii^* *^e..7:(,+,r)] • (77) 

Finally, we may come to the Hamiltonian. Let us assume for a moment that 
the loop we are dealing with has no intersection nor kinks, so that we can set 
the Einstein term to zero. Inserting our last result into the Hamiltonian we have 



F*[a]= lim Y^L^JjjS d3xJ2 S^iae,x)^;S -^[a] . (78) 

For L small enough, and assuming that S << L, so that "boundary effects" of 
the box can be neglected, every cube / contains only one end-point, and since 
Cwcyi gives zero unless there is an end-point in /, we have 



ma] = ^ hm J'-j^L^SY^ m i *H. (79) 



where we have restored the explicit expression for CrS- 

It is now time to study the limits explicitely. Let us first focuss on the crucial 
prefactor 

C{L,r,S) ^ l'^J-m=-^J^. (80) 

The question we have to address is the finitness of the above limit. Clearly the 
result depends on the order in which the limits are taken. This is precisely what 
we should have expected: different orders in which the limit is taken correspond 
to inequivalent definitions of the quantum operator. Since all these definitions 
correspond to the same classical limit, the choice amount to a choice of different 
orderings of the quantum hamiltonian. The question is whether there is one 
choice that gives us a finite quantum operator. 

Of course we may not confine ourselves to the choice between taking one 
first or another one first of the three limits: we choose any combination. More 
precisely, we may consider the three dimensional L,t,6 space, and study the 
limit of C(i, T, S) as we approach the point L = 0,t = 0,6 = 0: this point can be 
approached in a variety of alternative ways, not just along one of the coordinate 
axis. Let us introduce a parameter e, and consider a curve L{e), r(e), (5(e) in this 
tree dimensional space, such that £(0) = 0, r(0) = 0, S{0) ~ 0. Our problem is 
to understand wether we can choose this curve in such a way that the limit 

lhnC{L{e),T{e),S{e)) (81) 
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is finite. We are not free to choose the curve L(e),T(e), (5(e) in a completly 
arbitrary way, because there is a certain number of conditions that we have 
imposed on the regularization parameters along the way. First, of course, we 
must have L > 0, r > 0, (5 > 0. Then, have required r > (5, and, in order to 
avoid "boundary effects" in the box, L > S. Can aU the conditions be satisfied 
and a finite hmit be obtained ? 

The first crucial point to be noted is that powers of lengths cancell exactly 
and the quantity C'{L,t,6) is dimensionless. This is a necessary condition for 
having a finite limit. (It was precisely the fact that we got a divergent quantity 
with the dimensions of a length, measured in the background metric, that pre- 
vented us from achieving a background independent renormalization in section 
(5)). By itself, however, this fact does not suffice to garantee a well defined 
limit. We now claim that this limit can indeed be chosen consistently with all 
the demands, as follows 



Lie) = 


= ke^ 


r{e) = 


- ea. 


S{e) . 


. e^a. 



(82) 

where a is an arbitrary length, and fc is a arbitrary dimensionless positive num- 
ber. It is easy to see that all the requested conditions arc satisfied. In the limit, 
we have 



limC(L(e),r(e), J(6)) = . /^^ = -^ := A^ (83) 

£^0 2V7r y T(e)d^(e) 2^tx 

A is a free dimensionless constant that emerges from the regularization proce- 
dure. Thus, the prefactor is finite in the limit. 

Thus, we write the action of the hamiltonian as 

iJ*[a] = A^ XI (-^O "^ *["] (^^) 

where we have introduced the "end-point operator" 

.Fe*H = lim.Fj(''''(') *[a]. (85) 

We now examine this " end-point operator" , its action, its finiteness and its 
transformation properties under diffeomorphisms. 

Since b (the point splitting distance of the two grasps) goes to zero much 
faster than t (the length of the added loop, we can now simply take the (5 — > 
limit first, and the t ^ limit second. Let us consider the (5^0 limit of 
J"J''^[q;] (with finite t). It is easy to see that if the end- point is free, the action 
of the operator is simply to add a small straight line of length r = ea to the 
end point of the loop, in the direction of the loop tangent (goign out from the 
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intersection). If the end-point is not free, the action of the operator produces 
one term for each component of a emerging from the end-point. The terms 
corresponding to the component of a that ends in ae is again just an addition 
of a small straight line to the end-point, emerging from the end point in the 
same direction as the original loop; while the other terms imply the addition 
of the small loop and also a rerouting through the intersection. The rerouting 
pattern can be calculated in straight forward way from equation (|77|). 

Before taking the limit r ^ 0, let us now assume that ^[a] is a diffeomor- 
phism invariant state. Thus ^[a] depends on the diffeomorphism equivalence 
class of a only. If ae is a free end-point, we have then 

lim^:"vl/[a] ^ hm 2vl/[a * *Yc.^,^\J = 2* [a], (86) 

because for small enough r the added loop will not interject any other loop, and 
the addition of a small line at the end of a loop does not change the diffeomor- 
phism equivalence class of the loop. 

If, on the other way, ae is not a free end-point, then the loop a * *7'^ - does 
belong to a different knot class than a. For instance, assume the end point ae 
falls over a smooth component (3 of the loop a. One of the terms in (^) adds 
a small straight line to a. Since this in the direction in which a emerges from 
the intersection, its action is essentially to cut away an infinitesimal portion of 
the loop at the end point. This can be seen by noticing that we can assume 
without loosing generality that a is a straight line at the intersection, and using 
the retracing identity, to " retrace" a away from the intersection along the small 
added loop. In the limit, we go back to a; however, for every finite value of 
the regularization parameter, the knot class has changed: we are in the knot 
class in which a does not touch /3. Thus, by continuity in the limit we still 
have the knot class in which a does not touch f3. The key point, now, is that 
in any case, since \l/[a] is diffeomorphism invariant, for small enough t we have 
that J^j'^\E'[a] becomes independendent from r. Therefore the limit is the limit 
of a constant function, and therefore is finite. The action of the Hamiltonian, 
therefore, has opened up the intersection between (3 and the end point of a. 

It is clear that the resulting action of !Fe is well-defined on the diffeomorphism 
invariant states. Thus, the operator H is finite and diffeomorphism invariant in 
the limit. 

If we now reinstate CEintsoin y^ 0, we have 



H= Yl \jM, + XFe, (87) 



intersections i 
end— points e 



where M was constructed in Q. iJ is a finite operator defined on knot states. 
Its action follows immediately from the construction above. 



26 



5.3 Topological Feynman rules 

It is clear from the discussion above that the action of the Hamiltonian con- 
straint on the generahzed knot classes can be represented in a fully geometrical 
way in terms of the action of the operator on single intersections. The operator 
"opens up" intersections, and rearrange the roofings through the intersections. 



The details of this action can be computed from equations (773590), which we 
report here for easy reference. The finite geometrical action on the loops is 
defined by the operator ( |77| ) 

This operator attaches to every intersection and to every line coming out from 
an intersection two small straight loop, one for each one of the two possible 
rearranging of the roofings at the intersection between the added little loop and 
the line. Next, we take the limit ( p5| ) 

.Fe*[a] = lim.Fj(^)'^(^) *[«]. (89) 

And finally, the Hamiltonian is given by (BO) 



intersections i 
end — points e 

The matrix elements of the operators Mi ( js) ) and Fe can be directly computed 
between any two given knot states. The calculation amounts in a straightfor- 
ward exercize in geometry and combinatorics, starting from the two equations 
above. The next problem is to compute the square root of the resulting (infi- 
nite) matrix. We expect that the square root can be computed order by order 
as the complexity of the knots considered increase. Work is in progress to com- 
pute explicitely the matrix elements, and thus understand if the structure of the 
resulting matrix allows a simple argoritm for extracting the square root. The 
resulting geometrical action represents the equivalent of the Feynman rules for 
this diffeomorphism invariant, or (infinite dimensional) topological, theory. We 
suggest to denote them as topological Feynman rules. 

6 QGD: dynamics 

We are now in the position of describing the general structure of Quantum Grav- 
itational Dynamics, or QGD, the quantum theory of gravitationally interacting 
fermions, evolving in the clock time defined by a scalar field. 

A physical quantum state \K) of the theory is specified by a generalized 
knot, namely an open braid K of order N (with N open end-points, TV even), 
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with an arbitrary finite number / of intersections. A more accurate notation 
for these states is given in equation (^5|). The quantum dynamics is given by 
the matrix H in braid space, given in equation (po[), the matrix elements of 
which are computed, order by order, according to the geometrical and algebraic 
rules given by equations ( ^ and (|85|), and in ref.Q. We can interpret the 
matrix elements of H as first order transitions amplitudes in a time dependent 
perturbation expansion in the clock time T. In principle, the exponentiation of 
the H action gives the full evolution. 

6.1 The simplest states 

For instance, we can start from the simplest state formed by a single nonself 
intersecting open line. In terms of the notation (Eq), this can be denoted as 

|2,2,0;1°) (91) 

where we have indicated the simplest value of /C§, a single line, by Ij. The 
moduli space of free open-ends is clearly formed by a single point, and thus we 
do not need a.i parameters. 

There are two fermions in this quantum state. We have that Mi|2, 2, 0; 1°) = 
because there are no intersections of kinks in |2,2, 0; I2). On the other side, 
we have from ( pq ) that 

i^e |2,2,0;1^) = A2|2,2,0;1°), (92) 

so that we get 

7^|2,2,0;lO>=2A|2,2,0;l°). (93) 

Therefore |2,2,0;l2) is an eigenstate of the theory, or equivalently, the time 
dependent Schrodinger quantum state 



|2,2,0;1°, D^expjiA^^T^ |2,2,0;1°) (94) 

is a solution of the exact quantum interacting theory. (We have restored phys- 
ical units, for clarity.) Perhaps this state corresponds to an extremely simple 
"universe" in which there are only two fermions gravitating around each other 
in the simplest of the quantum geometries. It is suggestive to think at this state 
as a kind of " atomic" " ground state" of a simple 2-fermions universe. 

Next we can consider the generalized knot formed by n non intersecting 
copies of the above, and denote it as |2n, 2n, 0; 1%^- It is then straight forward 
to see that the time evolution of this state is given by 



2n,2n,0;l°„, T) = exp <^ z nA W— T I |2n, 2n, 0; 1°„) (95) 
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and that the corresponding energy eigenstates are 

En = nEi=nxJ^. (96) 

V Cj 

As soon as we consider simple intersecting states, the full complexity of the 
operators Mi and Fe becomes relevant, and we have non-trivial time evolution. 

6.2 Comments 

Before concluding, we list here a certain number of comments and considera- 
tions, as well as pointing out several important problems that remain open. 

1. Conservation of particle number. The operator F,, defined above acts on 
end-points by displacing them, and possibly by changing the associated 
roofings at intersections, but never creates or destroys end-points. Since 
the operator Mi too, conserves the number of end-points, it follows that 
the hamiltonian that we have defined conserves the number of particles. 
This is at first surprising, given that in general there is particle creation 
from space-time dynamics. But a more close analysis shows that this 
conservation is to be expected. Unlikely the Einstcin-Dirac theory, indeed, 
the Einstein- Weyl theory does conserves particle number. This can be seen 
classically from the fact that the quantity 



N := / d'^x V^(a;) nA{x) (97) 

commutes with all the constraints, including the hamiltonain constraint 



30 . One can immediately define the corresponding operator (say, using 



the Loop Transform), which turns out to be 

N\\N,I,D; ar....aT'; 1C°^ J = N\\N,I,D- a^.-ar; ^|^ ™> 

. . '^^^) 

This confirms the interpretation of the number N of end-points as the 

particle number. Since [N, H] = 0, the number of end-points iV is a 

conserved quantum number in the theory. 

Particle anti-particle distinction. The Weyl field theory describe a particle- 
antiparticle couple (say a neutrino and its anti- neutrino). In the Lan- 
grangian formulation the fermions are described by two complex fields. 
Since the action contains only first derivatives, the phase space has the 
same dimension as the space of the lagrangian fields, namely four real 
dimensions per point. These give two degrees of freedom, which describe, 
indeed, the particle and its antiparticle. Do the end-points of the loops 
represent particles or anti-particles ? The answer is that the distinction 
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is not gauge invariant, thus the question is not weU posed in the the- 
ory. In flat space one can globafly distinguish particles from antiparticles; 
but when the Weyl system is coupled with gravitation, something curi- 
ous happens: the particle antiparticle distinction becomes local. Consider 
two field excitations in two different space position, and assume the first 
is a particle; then, the fact that the second be a particle or an antipar- 
ticle depends on the parallel transport operator between the two. This 
is because the particle and the antiparticle are distinguished by different 
directions in the internal spin space, and we can only compare directions 
in spin space in different points by using the connection. Since the par- 
ticle anti-particle distinction is gauge dependent, in a gauge-independent 
formulation there is not way to distinguish particles from anti-particles. 
This is why the end-points of the loop represent at the same time both 
kinds of excitations. 

Regime of validity of the formalism, and complex energy eigenvalues. This 
is an important feature of the clock field formalism that we must be dis- 
cussed in detail. The formalism cannot be used in any regime of the 
system. This is already obvious at the classical level: Consider an arbi- 
trary solution of the field equations A{x,t),'ip{x,t),T{x,t): in general it 
is not possible to invert T{x,t) -^ t{x,T) globally. Thus, we certainly 
cannot use T as a time variable for every solutions of the field equations 
and for every spacetime region. On the other side, there are solutions 
and spacetime reions where we can make the inversion. Consider an ini- 
tial configuration of the A{x),tp{x),T{x) fields and their time derivatives 
A(x), "tpix), T{x) on a given space like surface. This defines a point in phase 
space. Assume that T{x) < in a spacial region R. There will be a time 
interval At for which T{x) will remain positive in R. More precisely, we 
can determine a region T in phase space, and a corresponding spacetime 
region S, such that for any initial condition in F, T{x) is positive in S. We 
shall say that the gravitational-fermion-scalar field system is in the clock 
regime in S if the initial conditions are in F, namely if T{x) is positive 
all over S. By definition, we can perform the inversion T{x, t) — > t{x, T) 
in the region in which the system is in the clock regime. Thus, as far as 
the classical theory is concerned, the formalism makes sense only in this 
regime. The same is true in the quantum theory. The quantum formalism 
that we have constructed is meanigfuU in the clock regime. 

A paradigm for this construction can be found in the quantum system of 
two uncoupled simple harmonic oscillator variables g{t), f{t), if we fix the 
total energy E and decide to never consider the evolution in the external 
clock time t, but rather use one of the two variables, say / as internal 
clock; namely if we decide to ask questions concerning what is the position 
g of the first oscillator, when the second is in /. We obtain the classical 
evolution g{f) by inverting f{t) -> t{f) and defining g{f) = g{t{f)). 
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We can also do the same in quantum mechanics (see ref. [Ell, where the 
example is worked out in detail). However, along any orbit there is a point 
in which the internal time variable / " comes back" (in t) , and therefore we 
obtain a non unitary evolution operator in /. The physical interpretation 
of this non-unitarity is clear: there is "no system" anymore for / arbitarry 
large. 

The formalism reflect this fact, both classically and quantum mechanically, 
in the form of the hamiltonian. The hamiltonian that evolves the two 
oscillators in the external time t is p^ + pi + g^ + ,p. The hamiltonian 
that evolves the systemin the internal time / is casilly obtained solving 
ioipf (see ref. Ill) 



Hif) = ^E-pl-g^-p (99) 

(where E is the total energy of the system in the t time). This (time 
dependent) hamiltonian generates the evolution equations for g{f). The 
important point to note is that the hamiltonian becomes immaginary for 
large /. This simply signal the fact that it doesn't make sense anymore to 
evolve g{f) in /. If we want to continue the evolution by using an internal 
time, we have to choose another, distinct, internal time, and "patch" the 
evolution. Note that this does not mean that the formalism that evolves 
in / is inconsistent or incorrect: it simply means that is has a certain 
domain of validity. The same holds in quantum mechanics. Indeed, it was 
shown in ref. E^ that the quantum hamiltonian corresponding to (M) is 
self- adjoint when suitably restricted to an (/—dependent) region of the 
Hilbert space, but develops immaginary eigenvalues if applied outside this 
region. 

Similarly, we expect that the the hamiltonian that we have defined will 
also have immaginary eigenvalues. These simply signal that twe are go- 
ing out from the domain of validity of the formalism, namely from the 
clock regime. We are trying to use evolution in T to describe the grav- 
itational field in regions where the T fields fails to be monotonic. Ex- 
plicitely, this possibility can be easilly traced back to the classical hamil- 
tonian constraint. Roughly speaking, since the form of this constraint is 
n^ -|- CEinstoin + C'weyi = (11 being the scalar field momentum) and since 
in order T to change sign 11 must vanish, it follows that the save region, 
is where CEinstcin + C'weyi > 0, which is of course a sufficient condition for 
the Hamiltonian we have definined, H = J A/CEinstein + C'weyi to be real. 
Thus, imaginary eigenvalues of H signal that we are exiting the regime 
of validity of the formalism we have developed here: the object we have 
chosen as clock is running backward. We must therefore exclude from the 
state space of QGD, as formulated here, the graphs that are eigenvalues 
of H witha an eigenvalue that is not a real positive number. 
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In particular, all the vacuum solutions of pure quantum GR that where 
previously found lie outside the clock regime. They are eigenstate of H 
whith vanishing eigenvalue. Classically, the vacuum solutions of the theory 
are only compatible with 11 = 0, namely T = 0, which clearly indicates 
that we are outside th eregime in which we can take the scalar field as a 
good clock. 

We do not consider this necessary restriction of the formalism as a serious 
limitation. Our long term aim is to develop a usable theory that can be 
employed, at least in principle, to describe Planck scale measurements 
and the Planck scale evolution of quantum geometry. We would be very 
content of having a sensible general covariant field theory that correctly 
describe this physics in the regime in which whatever we are using as a 
clock keeps behaving as a clock. 

4. Scalar product. One of the weak points of the Loop Representation is 
given by the fact that a complete and consistent definition of the scalar 
product is not yet available. The conventional wisdom is that once phys- 
ical observables on the physical state space have been constructed, the 
scalar product is determined by the requirement these physical observ- 
ables be self-adjoint. The present work is a step in this direction. The 
(real eigenvalues) states of the hamiltonian H must form, if the formal- 
ism is consistent, an orthogonal basis. Thus, working out explicitcly the 
eigenstates of H in knot space should at the same time lead to a partial 
definition of the scalar product. For a discussion of the definition of the 
scalar product for fermion in functional representations, see p5[ . 

5. Taking limits on knot space. Finally, let us discuss a subtle point in the 
definition of H , which we percieve as the most delicate and potentially 
problematic point in the construction above. We refer to the different 
way in which the (5—^0 and the r — > limit have been dealt with, when 
dealing with knot states. 

To focus the point, let us consider a model example. Consider the space 
C[B] of the continuous functions f{x) on the real line. Consider the closure 
D of the space C[B], say in the pointwise topology, such that D contains 
also piecewise continuous functions as the step function 9 defined by 

e{x) = 1 if X > ; e{x) = otherwise. (100) 

Now define the linear functional k on C[B] as follows: (fc, /) :— limx^of{x), 
and assume you want to extend k from C[R] to D. There are two possible 
strategies: one is to keep the definition 

(fc,/)=lim/(x). (101) 

x—>0 



32 



The other is to note that an equivalent definition of k on C[R] is (fc, /) = 
/(O), and thus to define 

{k,f)=f{0). (102) 

According to the first definition 

ik,0) = l, (103) 

according to the second 

{k,e) = 0. (104) 

We are in a similar situation when we need to study the action of the 
hamiltonian H on knot states. In quantum mechanics, operators are often 
defined on dense subspaces of the state space. For instance we begin by 
defining the momentum operator in Schrodinger mechanics not on the full 
Z/2 state space but on the dense subspacc of the differentiable functions; 
then we can extend it. The Hamiltonian H that we define in this paper 
contains a certain number of limiting procedures. We may first rigourously 
define it on a suitable restricition of the space of the loop functionals. For 
instance we may assume that the loop functionals are continuous in all 
the deformations that we consider. H is well defined on this space. Then, 
however, we want to consider the action of H on the knot states. These are 
not continuous in the deformations that we consider and thus we need to 
define a suitable extension of the operator. At this point we have a choice 
that essentially reflects the choice we described in the simple example 
above. As far as we understand, this choice, if not dictated by internal 
consistency, is again part of the quantization ambiguities as the ordering 
of the dynamical operators. 

The important point we want to make here is not that a choice of the 
extension has been made in computing the action of H on the knot states, 
but that two different choices have been made for the two limits (5^0 
and r ^ 0. In fact, as far as the (5 ^ is concerned, we have assumed 
that we should first take the limit, and then consider the extension of the 
action of the operator to diffeomorphism invariant states; while as far as 
the r — > limit is concerned, we have assumed that we should first extend 
the the operator to diffeomorphism invariant states, and then take the 
limit. 

This choice is not completely arbitrary: 5 must go to zero faster that r, 
and, if we take away the fake dimensions added by the integration, we see 
that the first significative term, which is the one that we are considering, 
is of order zero in 5 and of first order in r. Ths means that already at 
the classical level what we are doing is precisely considering a function 
/((5, r) and picking up terms of the form ^/((5, t)| 5=7-^0 • Thus it is not 
completely unreasonable that this difference gets translated in the differ- 
ent ways in which the two limits are taken on loop space: roughly, we are 
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"really" looking at the 5 — point, and we are "really" looking at the 
limit in the first order expansion in r. However, these are hand-waving 
justifications of our choice. Untill a well-defined calculus on Loop Space is 
constructed |3^, we do not see a way to transform these tentative explo- 
rations into solid mathematics. Our only real justification at this point, if 
any, is the hope that the (finite) structure we are constructing be internally 
consistent and, perhaps, related to Nature. 

6.3 What next 

The next step in the construction of the theory should be to compute explicitely 
matrix elements of Mj and Fg (see eq.(pO|)), and to understand whether there 
is a direct algoritm for extracting the square root. If this can be done, the 
theory is essentially at the stage where the evolution of physical states can be 
described. 

As noted above, the scalar product is partially fixed by the construction 
itself. The energy and the particle number are conserved observables. There 
are other observables in the theory that one may consider, and evolve, as the 
area observable discussed in references [|| Q . A crucial test for the consistency 
of the scheme developed here is, as was noted in [g|, whether the second order 
term of the time dependent perturbation expansion develops divergences. 
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